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ABSTRACT 

We report on the detection of variability from the proposed IR counterpart to the Anomalous X-ray Pulsar 
(AXP) IE 1048. l-5937based on Chandra and ESO opticaPIR deep observations earned out in 2001-2002. Within 
the narrow Chandra uncertainty region for IE 1048.1-5937 we found only one relatively faint (7=22. 1±0. 3, 7- 
Ks=2A) source, while the recently proposed IR counterpart was not detected down to a limiting Ks magnitude of 
^20.7 (3(T confidence level). This implies a remarkable IR brightening of this object, /S.Ks>\3, on a timescale 
of about 50 days. Although our knowledge of the IR properties of AXPs is rather limited (there is only another 
source, IE 2259+586, for which IR variability has been detected), the observed IR variability of the proposed 
counterpart strengthens its association with IE 1048.1-5937. Our results make the IR (and presumably optical) 
variability a likely common characteristic of AXPs, and provide new constraints on this class of objects. 

Subject headings: stars: pulsars: general - pulsar: individual: - IE 1048. 1-5937- infrared: stars - stars: 
variables: other - X-rays: stars 



1. INTRODUCTION 

After more than 20 years since the discovery of X-ray pul- 
sations from the prototype source IE 2259+586, the current ob- 
servational properties of AXPs are not yet conclusive to un- 
ambiguously assess their nature. There are currently five con- 
firmed members of the AXP class plus two candidates. Al- 
though we can be reasonably confident that AXPs are magnetic 
rotating neutron stars (NSs), their energy production mecha- 
nism is still uncertain. It is also unclear whether they are soli- 
tary objects (either low magnetised or magnetars) or are in bi- 
nary systems with very low mass companions (for a review see 
Israel et al. 2002; Mereghetti et al 2002 and references therein). 
Different production mechanisms for the observed X-ray emis- 
sion have been proposed, involving either accretion or the dissi- 
pation of magnetic energy. The recent detection of X-ray bursts 
from IE 2259+586 and IE 1048.1-5937 has strengthened the 
possible connection of AXPs with Soft Gamma-ray Repeaters 
(SGRs; Kaspi & Gavriil 2002; Gavriil et al. 2002). More- 
over the detection of relatively large pulsed fraction pulsations 
in the optical flux of 4U 0142+614 (at the same period of X- 
ray pulses) seems to rule out the X-ray reprocessing and fa- 
vors models based on isolated NSs (Kern & Martin 2002). The 
magnetar model, originally proposed by Thompson & Duncan 
(1993) to explain the properties of SGRs, appears currently to 



be fairly succesful at interpreting these properties of AXPs. 

As done in the past in the case of High and Low Mass X-ray 
binary system (HMXBs and LMXBs), finding the optical/IR 
counterparts to AXPs might give the key to definitively assess 
the nature and the emission mechanisms of AXPs. The lack of 
any known distinctive property of AXPs in the optical/IR bands, 
makes the X-ray positional accuracy at the (sub)arcsec level a 
fundamental preliminary step to carry out very deep searches 
for potential opticaPIR counterparts. Moreover a high X-ray 
spatial resolution can be exploited to search for small scale dif- 
fuse emission or the presence of structures in the vicinity of the 
sources. 

IE 1048.1-5937 has been extensively studied in the X-ray 
band since its discovery in 1979 (Seward et al. 1986). Previous 
searches for an optical counterpart gave negative results (Se- 
ward et al. 1986; Mereghetti et al. 1992; Wang & Chakrabarty 
2002a). Recently, a relatively faint IR counterpart (7=21.7; 7- 
Ks=2.y) to IE 1048.1-5937, with unusual IR colors, has been 
proposed (Wang & Chaki-abaity 2002b; hereafter WC02). 

In this paper we report on the early results obtained for 
IE 1048. 1-5937 as a part of a joint ESO/Chandra large project 
aimed at the identification and study of the optical/IR counter- 
parts of AXPs in the southern sky, and the study of the spa- 
tial distribution of X-ray emission from AXPs (results con- 
cerning Chandra observations of IRXS J 170849-4009 10 and 
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AX J 1844-0258 will be reported elsewhere). 

2. Chandra OBSERVATIONS 

The field including the X-ray position of IE 1048.1-5937 
was observed during Chandra Cycle 2. Previously determined 
positions of the source had limited accuracy; 15" for EXOSAT 
(Mereghetti et al. 1992) and ^4" for XMM-Newton (Tiengo 
et al. 2002). The source was observed with the High Reso- 
lution Imager (HRC-I; Zombeck et al. 1995) in order to ob- 
tain the most accurate X-ray position. IE 1048.1-5937 was 
observed with Chandra on 2001 July 4 for an effective expo- 
sure time of 9870 s. Data were reduced with CIAO version 2.2 
and analysed with standard software packages for X-ray data 
(CIAO, XIMAGE, XRONOS, etc.). The observation was car- 
ried out with a nominal aspect solution and latest calibration 
files were used. Four different detection algorithms were used: 
the CIAO wavdetect and celldetect tool, the XIMAGE 
sliding cell, and the wavelet-based pwdetect software de- 
veloped and optimised for Chandra images (Damiani et al. 
1997). Based on a statistical study performed on several point- 
like sources detected within the field of view of public HRC-I 
archival data with known optical counterparts, we found a dif- 
ference up to 0."2 in coordinates among the four different de- 
tection methods. In the following, we refer to the coordinates 
obtained with the latter method (we assumed that differences in 
the coordinates among the four methods are purely statistical). 
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Fig. 1. — Radial surface brightness profile of the Chandra HRC-I emis- 
sion from a circular region of 7" radius around IE 1048.1—5937. The dashed 
line represents the simulated Chandra point spread function for a point-like 
on-axis source. 

Only one source was detected in the HRC-I. This has 
the following coordinates: R.A. = 10*'50"07?12, DEC. = - 
59°53'21."37 (equinox J2000), with an uncertainty circle ra- 
dius of 0"7 (90% confidence level). Note that the Chandra po- 
sition is not consistent with that obtained with XMM-Newton 
(Tiengo et al. 2001; R.A. = 10*'50'"06.^3, DEC. = -59°53'17"; 
equinox J2000, 90% confidence level circle radius of 4"). How- 
ever, inspection of the attitude data during the XMM-Newton 
observation, showed that an error radius of about 9" is likely 
more realistic (Tiengo 2002, private communication). Photon 



arrival times were extracted from a circular region centered on 
the previous reported position with a radius of 1''5, including 
more than 90% of the source photons. These were corrected 
to the barycenter of the solar system and timing analysis ap- 
plied (CIAO f xbary tool). A coherent pulsation at a period 
of 6.4529±0.0001 s was detected confirming that the detected 
source is indeed the AXP IE 1048.1-5937. The latter value 
falls within the expected period interval extrapolated from the 
more accurate timing analysis obtained from a long-term /?XTE 
monitoring program of IE 1048.1-5937 (Kaspi et al. 2001; we 
assumed the maximum and minimum P reported in their Ta- 
bid). 

One of the aims of the Chandra HRC-I observations of 
AXPs, was to study the radial profile of the X-ray emission 
from the AXPs with the unprecedented spatial resolution of- 
fered by the HRC-I. Figure 1 shows the spatial profile we ob- 
tained for the (background subtracted) X-ray emission from 
a region of 7" radius centered on the above reported source 
position. The radial profile is in good agreement with the ex- 
pected Chandra point spread function (solid line simulated with 
MARX; Wise et al. 1997). 

3. ESQ OBSERVATIONS 

Over the last two years we carried out deep optical and IR 
observations of the field including the original Einstein po- 
sition of IE 1048.1-5937. These data were obtained from 
the 3.5-m New Technology Telescope (NTT, La Silla, Chile) 
equipped with the SUperb Seeing Imager-2 (SUSI2; 2 
CCDs of 2048x4096 pixels; pixel scale of 0"16 and 5.5' x 5.5' 
field of view) in the optical band, and with Son OF Isaac 
(SOFI; Hawaii HgCdTe 1024x1024 aiTay; 0."292 pixel scale 
and 4.9' x4. 9' field of view) in the near-IR band. 

Optical Johnson filter deep images (BVRI) were first car- 
ried out with BUS 1 2 in service mode on 2001 March 20 with 
the following effective exposure times: 2700 s (B), 4350 s (V), 
2400 s (R) and 3400 s (/). Weather conditions were good with 
seeing of about 1". Standard reduction packages were ap- 
plied to the optical data (DAOPHOT II; Stetson 1987) to obtain 
the photometry of each stellar object in the images. Limiting 
magnitudes of 25.5 (B), 25.5 (V), 24.5 (R) and 23.5 (/) were 
reached. 

JH images were initially carried out with SOFI in service 
mode on 2001 May 29; 2700 s (each filter) of effective exposure 
time. Unfortunately, during observations clouds were present 
and the seeing was poor (1"-1."5) yielding a limiting magni- 
tude (3ct confidence level) of 22.5 (J) and 21 (//). SOFI JHKs 
images were obtained with better weather (photometric) and 
seeing (0"6) conditions on 2002 February 19. Final exposure 
times were of 3600 s (J), 4200 s (H) and 4800 s (Ks). Single 
10 s-long exposure images were taken for each filter (5 s-long in 
Ks) with offsets of 40" in order to sample and subtract the vari- 
able IR background. Standard IR software packages were used 
for sky frame subtraction and image coaddition (Eclipse and 
IRDR; Devillard 1997 and S abbey et al. 2001). During the 
latter run limiting magnitudes of 7^23, i/^21.5, and Ks^lQ.l 
were reached. 

To register the Chandra coordinates of IE 1048.1-5937 on 
our optical/IR images, we computed the image astrometry us- 
ing as reference the positions of stars selected from the USNO- 
A2.0 catalog which has an intrinsic absolute accuracy better 
than 0."25 (depending on magnitude and sky position of refer- 
ence stars; Assafin et al. 2001). About 40 objects with mag- 
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Fig. 2. — Optical / band (left panel) and near-IR Ks band (right panel) images of the region including the position of IE 1048.1—5937, taken from the 3.5-m ESO 
NTT. The Chandra uncertainty circles are also shown: the thick solid line circle has radius 0" 8 and represents the position reported in this paper, while the thin one 
is the WC02 position with a radius of 1". The faint object marked with X3, which lies on the upper edge of the circles, is the only one consistent with the X-ray 
position. The square represents the position expected of XI, proposed as the IR counterpart to IE 1048.1-5937 by WC02. Coordinates are RA (hms) and Dec. (" ' 
"; equinox J2000). 



nitudes brighter than R=14 were identified in the optical and 
IR images and used for astrometric calibrators. The pixel co- 
ordinates of the reference stars were computed by a two di- 
mensional Gaussian fitting procedure and transformation from 
pixel to sky coordinates was then computed using the Graphical 
Astronomy and Image Analysis Tool (GAIA."). After taking 
into account the uncertainties of the source X-ray coordinates 
(Of'7), the rms error of our astrometry (0."15) and the propaga- 
tion of the intrinsic absolute uncertainties on the USNO-A2.0 
coordinates (0"25), we estimated the final accuracy to be at- 
tached to the IE 1048.1-5937 position as about 0."8. Figure 2 
shows a region of 20"x20" ai'ound the IE 1048.1-5937 posi- 
tion in the / and Ks bands with the Chandra uncertainty circle 
superimposed (90% confidence level). 

4. DISCUSSION 

The NTT IR data show that no object is consistent with the 
Chandra position of IE 1048.1-5937 in the optical band down 
to a y limiting magnitude of about 25.5 (see Figure 2). On the 
other hand, a faint object (marked X3 in Figure 2) is present 
in the IR Ks frames (7=22.1±0.3, Ks=\9.1±Q3). This object 
is also detected in the WC02 frames, although its magnitudes 
are not listed. A visual inspection of the WC02 image sug- 
gests that the source X3 has remained nearly unchanged in flux 
with respect to our frames (superimposing to our IR images the 
Chandra position and uncertainty radius reported by WC02 for 
IE 1048.1-5937, the object X3 faUs within the circle). Note 
that during the NTT IR observations carried out on 2001, the 
object X3 is not detected due to the poor seeing conditions and 
contamination with the nearby bright object. More surprisingly, 
we note that no object is detected at the position of the source 
XI reported by WC02 (the position of XI is marked with a 
square in Figure 2). For comparison, the object marked X2 in 

" See http://star-www.dur.ac. uk/'^pdraper/gaia/gaia.html 



Figure 2 has a magnitude of y=20.63±0.06, //=19.73±0. 15 and 
A'i=19.0±0.2 (consistent, within 1-1.5(7, with the values re- 
ported by WC02). As the limiting magnitudes of our IR frames 
are well above the values measured by WC02 for XI, we con- 
clude that the source was considerably fainter during our ob- 
servations. This IR brightening {/S.Ks >1.3; about a factor of 
3 in flux), displayed by the source on a timescale of about 50 
days (the time interval between our last IR images and those of 
WC02) further strengthens the link between the IR source and 
IE 1048. 1-5937. Moreover, the chance probability of finding a 
highly variable IR source within a circle of 0."8 is negligible. 

An IR counterpart was recently found for the prototype of the 
AXP class, IE 2259+586, which also showed IR variability at a 
levelof A/:.s--1.3(Kaspietal. 2002). This IR variability might 
be associated to the short-term X-ray bursts (similar to those 
displayed by SGRs) observed from the source three days before 
the IR observations (Kaspi & Gavriil 2002). Recently, short X- 
ray bursts were detected in the X-ray flux of IE 1048.1-5937 
(Gavriil et al. 2002). All these findings together, suggest a 
close similarity of IE 1048. 1-5937 and IE 2259+586, such that 
the IR counterpart of the former might be expected to be vari- 
able. For this reason XI is the most likely IR counterpart 
to IE 1048.1-5937, even though source X3 is included in the 
Chandra uncertainty circle and its color {J-Ks=2A) is within 
the value interval predicted for AXPs (see also Figure 2 of 
WC02). 

The fact that IE 1048. 1-5937 showed two X-ray bursts close 
to each other, implies that many other bursts might have been 
missed due to the sparse X-ray monitoring with the /?XTE. 
Moreover since our non-detection of the IR candidate predates 
the WC02 detection, it would appear more natural to asso- 
ciate (if at all) the IR variability to some X-ray activity that 
took place between the two IR campaigns (as opposed to the 
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few second-duration X-ray bursts that was detected 4 months 
earlier). Finally we note that the X-ray/IR observations of 
IE 2259+586 are such that we do not know whether the de- 
tected IR brightening follows or precedes that in the X-ray 
band. Future multi-wavelength monitoring programs will prob- 
ably allow us to better understand these points. 

The non-detection of XI in our IR frames implies that the 
"quiescent" (or low state) IR emission level of IE 1048.1-5937 
is still undetected. Correspondingly, it cannot be ruled out yet 
that the spectral flattening proposed by WC02 might result from 
wavelength-dependent variability. Similarly, the measured IR 
fluxes for the counterpart cannot be used to study the overall 
(from X-ray to IR) energy spectrum of IE 1048.1-5937, since 
they reflect a different energy distribution than the quiescent (or 
low state) one. Deeper IR observations are needed to determine 
the quiescent emission state of this source. 

IR variability might represents a new, fairly common, char- 
acteristic of AXPs (and possibly SGRs), and potentially a di- 
agniostic of their nature. We note that, so far, very little is 
known on the expected opticaWR emission from a magnetar. 
Both models based on magnetar and accretion from "fossil" 
disk have difficulties in accounting for transient and/or variable 
emission at wavelengths shorter than those in the X-ray band 
(Thompson & Duncan 1996 and references therein; Chatterjee 
& Hemquist 2000 and references therein). On the other hand, 
the analogy with the X-ray bursts of LMXBs, suggests that the 



effects of the X-ray burst of AXPs might propagate towards 
longer wavelengths (e.g., via reprocessing). It is somewhat in- 
triguing that at least the values of AKs and colors of the IR 
counterpart of IE 2259-H586 and IE 1048.1-5937 are similar to 
those displayed by some classes of accreting cataclysmic vari- 
ables and low mass X-ray binaries hosting a neutron star. 

The superb Chandra HRC-I spatial resolution allowed us to 
look for extended X-ray emission around the pulsar; this might 
result from an X-ray nebula powered by the AXPs as well as 
a scattering halo. A first detailed study was carried out for 
IE 2259+586, where emission beyond 4" and extending up to 
100" was found (Patel et al 2001). However in this case the 
source is embedded in the diffuse emission from the supernova 
remnant G109. 1-1.0, and it is very difficult to disentangle this 
component from that (possibly present) around the AXP. The 
present Chandra HRC-I data are by far the best spatial reso- 
lution observation carried out for an AXP and no significant 
diffuse X-ray emission was found beyond -^2". 
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